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Abstract. We discuss a very general theory of gravity, of which Lagrangian is an arbitrary function of the curvature invariants, 
on the brane. In general, the formulation of the junction conditions (except for Euler characteristics such as Gauss-Bonnet 
term) leads to the powers of the delta function and requires regularization. We suggest the way to avoid such a problem by 
imposing the metric and its first derivative to be regular at the brane, the second derivative to have a kink, the third derivative of 
the metric to have a step function discontinuity, and no sooner as the fourth derivative of the metric to give the delta function 
contribution to the field equations. Alternatively, we discuss the reduction of the fourth-order gravity to the second order 
theory by introducing extra scalar and tensor fields: the scalaron and the tensoron. In order to obtain junction conditions we 
apply two methods: the application of the Gauss-Codazzi formalism and the application of the generalized Gibbons-Hawking 
boundary terms which are appended to the appropriate actions. In the most general case we derive junction conditions without 
assuming the continuity of the scalaron and the tensoron on the brane. The derived junction conditions can serve studying the 
cosmological implications of the higher-order brane gravity models. 
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1. INTRODUCTION 

In recent years there have been growing interest both in brane universes HJEJ and in the higher-order gravity theories 
of which the simplest is f(R) gravity [3 1. In this talk we are going to combine both ideas and formulate the higher-order 
gravities on the brane. It emerges that the formulation is a bit non-trivial, since one faces ambiguities of the quadratic 
delta function contributions to the field equations. We will say how to avoid these problems and show how the Israel 
junction conditions for such higher-order brane gravity models can be formulated. 



2. FOURTH-ORDER GRAVITIES. 

When one considers the general gravity theories (e.g. flUEl): 

S = X~ l Jd D V=gf(X,Y,Z) (1) 
in a D-dimensional spacetime (x = const.), where X,Y,Z are curvature invariants 

K, I — K ab K , Z — K ahcd K , {Z) 

then one immediately faces the 4th order field equations, except when they reduce to the theories with Euler densities 
of the n-th order /W flg] 

S= ( Av^^^"', (3) 

JM „ 

the lowest of them being the cosmological constant I^°> = 1 (K"o = — 2A(2rc 2 )~ 1 = — 2A/167TG), the Ricci scalar 
/W = R (ffi = (2k- 2 )-'), and the Gauss-Bonnet density 1^ = RGB = R 2 - 4R ab R ab + R abcd R abcd (k 2 = a{2K 2 )-\ 
a = const.). 

However, the theories based on the Lagrangians which are the functions of the Euler densities such as 

f(R)=f(X), f(R GB )=f(Z-4Y+X 2 ) , /-/(/ (n) ) (4) 

are again fourth-order. 



3. FORMULATION OF THE 4TH ORDER GRAVITIES ON THE BRANE - ISRAEL 

FORMALISM. 



In the context of the recent interest in string/M-theory, it is interesting to formulate the general gravity theories ([T| 
within the framework of the brane models [7]. The full brane action for such a theory reads as 

S = X' 1 I d D xy/=if(X,Y,Z)+S bmne + S m , (5) 



M 

with the total energy-momentum tensor 

Ta = T a b -9(-w) + T a b +0(w) + 8(w)S a b , (6) 

where S b is the energy-momentum tensor on the brane, and T a b ± are the energy-momentum tensors on the both sides 
of the brane, 9(w) is the Heaviside step function, and 8(w) is the Dirac delta function. 
We assume Gaussian normal coordinates, i.e., (/i, V = 0, 1,2, . . . ,D — 2; w = D) 

ds 2 = g a i,dx a dx b = edw 2 + hfi V dx^dx v , (7) 

where £ = n ■ n = +1 for a spacelike hypersurface, e = — 1 for a timelike hypersurface, and h a j, = g a b — £n a rib is 
a projection tensor onto a (D — 1) -dimensional hypersurface, n is the normal vector to the hypersurface. In these 
coordinates the extrinsic curvature is 

1 dh^v 

K^iv — ^ ]5 — , (8) 

2 dw 



and the Gauss-Codazzi equations read [ 8 1 



Rwfiwv — ~~ 5 h KpvK fi, (9) 



dw 

Rwfivp — V v K il p~'VpK jl v, (10) 
R Xfivp = ^^RXfivp + e [K^vKxp -K^pKxv] ■ (11) 

In the standard Israel approach |9 1 one assumes that at the brane position w = 0: 

hjtv = h+ v , (12) 

hfiv,w 7^ n fiv,w > ^juv ^ ^fiv ' (13) 

i.e., the metric is continuous but it has a kink, its first derivative has a step function discontinuity, and its second 
derivative gives the delta function contribution. 

In terms of 9(w) and 8(w) functions this is equivalent to 

VvW = h- v (w)e(-w)+h+ v {w)6{w) , (14) 
dh „ v <5/i+ , . an uv , , 

= ^ie(-w) + ^d(w), (15) 

d 2 hp, v 
dw 2 

I dh7,„ dhf,,, \ 

(16) 

For the standard brane models with the Einstein-Hilbert action in the bulk 

S = / d D xy/^gR + S hrane + S m (17) 
2jc Jm 




the field equations read as O 



Gw 
w 



-\ {D - l) R+\e[K 2 -Tr(Ki)]=K 2 TZ, 

e[y n K-y v Ki] = K 2 r; 1 



dK% _ „ dK 

~~ 5 ^y- 

ow dw 



(18) 
(19) 

(20) 



-KK% + ^ v Tr(K 2 ) + ^ V K 2 



2 T H 



K T 



and in the limit lim H ^o J- w , which "fishes out" the delta function contributions, one gets the standard Israel junction 
conditions as (8) : 

e{[K^]-8^ v [K]} = kV v , 1JC%]=K%+-K%-. (21) 

By [X] = X + +X~ we define a jump of an appropriate quantity X at the brane. 

However, for the general f(X,Y,Z) theory on the brane, the standard continuity relations ( 12 1-( 13 i do not work. 
This can be seen from the field equations of the action ([T]) 



Pab = 



;T a b, 



(22) 



006 



- \fg° b + fxR ab + 2f Y R^"R b \ + 2f z R ed ^ a R b \. de 
fx-,cd{g ab g cd ~g ac g bd ) + n{f Y R" b )+g a \f Y R"% cd 
2{fYR c{a ) b) c -A(fzR d{ab)c ^ 



);cd, 



(23) 



where fx — d f/dX etc. 

Take, for example, the square of the Ricci scalar 



R 



D-l 



R + e 



■3Tr(K 2 )-K 2 



where K = , Tr(K 2 ) = K^K^y, and appropriately, of the Ricci tensor, and of the Riemann tensor. These squares 
produce the terms of the type 



d 2 h^ V dKf.y dKftydK^ fdK\ 



d 2 w dw ' dw d 



dw J 



(24) 



which are proportional to 8 2 (w), and so they are ambiguous. 

Amazingly, all these ambiguous terms cancel each other exactly in the case of the Euler densities [10]. In fact, the 
junction conditions for one of the Euler densities - the Gauss-Bonnet density, were already obtained as lfTTll2Tll 



where 



2a (3[7 MV ] - [J}h^ v -2[P} tipva [K}P°) + [K^] - [K]h^ v = -k 2 S^ 



Pfipva — R^ipva + 2/!jj[cj/?v]p + 2/Zp[ v /? CT j^ + Rh^h^p , 



(25) 

(26) 
(27) 



In the limit a — > 0, they just give Einstein-Hilbert action junction conditions (21 



In view of the ambiguities of the terms in ( |24) , we find two ways to formulate the junction conditions for general 
f(X,Y,Z) theories on the brane. 



A. Smoothing out the continuity conditions for the metric tensor at the brane 



In order to do that we impose more regularity onto the metric tensor at the brane position, i.e., we consider a singular 
hypersurface of the order three [9| which fulfills the conditions (compare (p"2|)-([T3]>) 
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A /IV,W "^iV.W 1 


fJ.V,WWW 
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(28) 
(29) 
(30) 
(31) 



i.e., the metric and its first derivative are regular, the second derivative of the metric is continuous, but possesses a 
kink, the third derivative of the metric has a step function discontinuity, and no sooner than the fourth derivative of the 
metric on the brane produces the delta function contribution. 

The physical interpretation as put in terms of the second-order theory can be that there is a jump of the first derivative 
of the energy-momentum tensor (e.g. jump of a pressure gradient) at the brane. 

In his seminal work, Israel [9| proposed a singular hypersurface of order two, which physically corresponded to a 
boundary surface characterized by a jump of the energy-momentum tensor (e.g. a boundary surface separating a star 
from the surrounding vacuum) which was characterized by 
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"/JV,w i 
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(32) 
(33) 
(34) 



i.e., the metric is regular, theirs? derivative of the metric is continuous, but possesses a kink, the second derivative 
of the metric has a step function discontinuity, and the third derivative of the metric on the brane produces the delta 
function contribution. 



The appropriate junction conditions can be obtained as follows. We rewrite the field equations (22 1-(23 i as 

X ^ab/ 



-gc ah w abd , d - 



C a bV- 



(35) 



where we have introduced is an arbitrary tensor field C ab , and 



W 



abd 



V 



ah 



fxAg ab g cd -g {ac g b)d )+(fYR ab y d 



g " b (f Y R^). c -2(f Y R^ a y b >-4(f z R d(al, > c ), c 



•ycd 



d{a\;b) 



ryd(ab)c\ 



i 



fg ab +fxR aa + 2f Y R c[a R 



c(a nb) 



2fzR edc(a R b) 



cde 



(36) 



(37) 



contain third derivatives of the metric giving a step function discontinuity, so that W -d is proportional to 8(w). 
Then, we integrate both sides of the formula (35i over the volume V which contains the following parts (cf. Fig. [TJ: 
G\,G1- are the left-hand-side and the right-hand-side bulk volumes which are separated by the brane, Al = dGl +A0, 
A2 = dG2 — AO are the boundaries of these volumes, and AO is the brane which orientation is given by the direction 
of the normal vector n. 
We have 



G1+G2 



- g c ab w abd . 4 da- 



G1+G2 



gC ab V ab d£l- 



^T lb c ab ^-gda , 

G1+G2 2 



(38) 



and so 



-g(c ab w abd ). d da 



G1+G2 

X rj,ab 
Gl 2 



-gdCl- 



G1+G2 

X rj^Ctb/^ 

G2 2 



gc ab . 4 w abd tm- 



-gdCl- 



X aab/^ 
TT^ C ab 
AO 2 



G1+G2 

yda , 



-gc ab v ab da 



(39) 



V 







G2 


G1 











FIGURE 1. A schematic picture illustrating the domains of integration used in derivation of the junction conditions. Here 
V = Gl + G2 is the total volume, Gl, G2 - are the left-hand-side and the right-hand-side bulk volumes which are separated by 
the brane, Al = dGl +A0, A2 = 9G2 — A0 are the boundaries of these volumes, and AO is the brane which orientation is given by 
the direction of the normal vector n. 



of which the first term can be integrated out to a boundary A1+A2 and then the limit V — > AO (or lim w _>o J- w m 
Gaussian coordinates) is taken. 

The final form of the junction conditions which generalize pT| onto the fourth-order gravity are 

[W] abd n d -*S ab = 0, [w] abd = W abd+ -W abd -. (40) 

It is remarkable that these junction conditions involve the higher derivatives of the scale factor. To see this take for 
example f(X,Y,Z) = f(R) theory in D = 5 dimensions with metric 

ds 1 = -dt 2 + a 2 (t,w)[dr 2 + r 2 (d® 2 + sin 2 ®d(l) 2 )} + dw 2 . 



The junction conditions (40 » give a jump of the third derivative of a(t,w), as expected 

[a'"] = |fl po, (41) 
Po = Po, (42) 
where (...)' = d/dw, do = a(w = 0), and the brane energy-momentum tensor is = (—po,po,po,Po)- 



B. Reduction to an equivalent 2nd order theory 

Yet another way to obtain the junction conditions is the reduction of the action ([T]i to a second-order action. This 
gives equivalent junction conditions, though at the expense of introducing a new tensor field H abcd (tensoron). In fact, 
starting from the action [ T2l 

Sg = X' 1 I d D x^gf( 8ab ,R abcd ). (43) 



IM 

we may transform to an equivalent 2nd order action in the form 

hi 



Si = / d u x^g{H^{R ghlj - f gUJ ) + f{ gab ^ cdef )} , (44) 



where 

d 2 f(gab,^abcd) 



H ghij _ df(gqb,§abcd) 



d(j) g hijd<j>klmn 



^0. (45) 



This transition for f(R) theory requires a new scalar H = f'(Q) (a scalaron) with the condition that f"(Q) ^ 0, and the 
equation of motion Q = R. Similarly, for /(Rgb) theory, one defines a scalar H — f'(A), with the equation of motion 
A = Rgb)- In order to get junction conditions, we have to slightly redefine the tensoron 

j^abcd }_^^acdb _|_ j_^abdc j^cbda jjacbd jjabcd _|_ j_^cbad^ (46) 

which in a particular case of f(X,Y,Z) theory takes the form 

piibcd _ uadgcb _ g cd g ha ) + f Y (2R ad g bc - R cd g ba - R ba g cd ) + 4 f z R acbd . (47 ) 

The field equations for an equivalent action d43) read as 



„ dV( gab ,H"*f) 

K shij Jjjghfj . (48) 

1 ab f + *L+H l » ed F ecd (g ab ,H kbm ) + {A^ d }, dc = -\r' b , (49) 

where 

V(gab,H cde f) = ^H h ^^ hij {g ab ,H cde f)+f{g ab , hlmn{gab,H cdef )) (50) 



In fact, the possibility to express the fields fygbcd as a function of g ab and H cde f is guaranteed by the condition ( 45 1 
(which is an analogue of the condition f"(Q) ^ 0). 

One can show that junction conditions of the second-order theory are equivalent to junction conditions of the fourth- 
order theory (7) . 

Applying the same method as in the previous case (i.e. taking the limit of V — > AO) we notice that the first three 
terms of ( [49] ) do not give any contribution to the junction conditions (since they do not contain delta functions at all) 
which now have the form: 

[A lab)cd ;d]n c = -\s" b . (51) 

Assuming that 

f{gabAabcd) = f^ab ab Aacb C r Ch c tabcd***) , (52) 

we can get the same result as in the 4th theory 

[A [ab)cd . 4 ]n c = [A^ cd . c ]n d = [-{fxAg ah S cd - g c(a g b)d ) 

+ (f Y R ab y d + g ah (fYR cd )-c (53) 
_ 2{fyR d(ayb) _ A{f z R d ^ c ), c }]n d = -[W abd ]n d . 

Similar approach was used for less-general f(R) theories of gravity on the brane by Borzeszkowski and Frolov 
El; Parry at al. OH, Deruelle et al. (BJ, and for f(X,Y,Z) = aX 2 + bY + cZ (a,b,c = const.) theories by Nojiri and 
Odintsov J5J. 



4. FORMULATION OF THE 4TH ORDER GRAVITIES ON THE BRANE - 
GIBBONS-HAWKING BOUNDARY TERMS 



In this approach, following the idea of Gibbons and Hawking ifTBI . we do not assume any vanishing of the first 
derivative of the variation of the metric tensor 8g ab ^ c on the boundary of the integration volume while using the 
variational principle. Strictly speaking, only the assumption of the vanishing of the normal derivative of the variation 
of the metric tensor 8g abw is required. Instead, we postulate that some extra terms to the action are added and that these 
terms "kill" the first derivatives of the metric variation. These terms are called Gibbons-Hawking boundary terms now. 
In fact, the Gibbons-Hawking boundary term for the Einstein-Hilbert action is composed of the trace of the extrinsic 
curvature and it was found by Gibbons and Hawking themselves (16). Then, for the action being the combination of 
the square of the Weyl tensor and an arbitrary function of the scalar curvature they were found by Hawking and Lutrell 
ifTTl and Barrow and Madsen [ 18"]. For the Gauss-Bonnet and other Lovelock densities they were found by Bunch 
[ 19 1, Mueller-Hoissen and Myers EOl . Davis lETl and Gravanis and Willinson [22|. The boundary terms for the action 
being an arbitrary function of the curvature invariants were found by Barvinsky and Solodukhin [23 1. 

For the theories which are of interest for this talk, the Gibbons-Hawking boundary terms have the following form 
11241 : for the f(R) theory the term reads as 

S GHp = -2{-l) p e [ yf^hHKd D - l x, (54) 

JdM p 

where H — f'(Q) is the scalaron, while for the f(X,Y,Z) theory it reads as 

Sgh, p = -(-1) p [ d D - l x^hA^ cd n c n d ^ ngab , (55) 

JdMp 

where A^ ab ^ cd is the tensoron. 

Using the method of the boundary terms we derived the most general Israel junction conditions for f(R) theory as 
USD: 

[K] = 0, (56) 
S ub n a n b = 0, (57) 
S" h h ac n b = 0, (58) 

-(D-l)[H :c n c ]-D[H]K = e^S ah h ah , (59) 

-h ab [H, c n c ]-[H]Kh ah + [HK ab ] (60) 

= e^S cd h ca h db . 

A generality of these conditions refers to the fact that no assumption about the continuity of the scalaron on the brane 
has been made. They reduce to the conditions already obtained in the literature, if one assumes [H] = [ 15 1. 

On the other hand, the most general Israel junction conditions for the f(X,Y,Z) theory, with no assumption about 
the continuity of the tensoron on the brane, are ll24l : 

[KA^ cd ]n c n d + [^ n A (ab ^ d }n c n d (61) 
- e [A^ d K cd ] - g" h [A^ d K ef ]n c n d 

+ 2e[D s A {ef ^ d n c n d ]h s e h ( "n b) - 2e[A^ cd . {c ]n d) = |s flfe , 

"6"c [^agcul] - n u n c [3(ngdb] - n b n d [^ngac] + n a n d [^ngcb] = . (62) 
They reduce to the conditions ( [51} , if one assumes continuity of the tensoron on the brane 

[A^ cd ] = . (63) 



5. FOURTH-ORDER GRAVITIES AND STATEFINDERS 



We claim the fact that general f(R,R a i,R ab 7 R a bcdR abcd ) theories are fourth-order may have some advantageous conse- 
quences onto their observational verification by the application of statefinder diagnosis of the universe. 

In fact, statefinders are the higher-order characteristics of the universe expansion which go beyond the Hubble 
parameter H and the deceleration parameter q: 



a 1 a da 

H = — , 9 = --2_ = --2 . (64) 
a H L a a L 



They can generally be expressed as (i > 2) 

* = <-^=<-.)^, 

and the lowest order of them are known as: jerk, snap ("kerk"), crack ("lerk") 

1 d da 2 1 d dc? 1 a^a 4 

j = fPa = ^' lT 4 ~a = ~lF ' fP^ = ^^ 1 ( } 

and pop ("merk"), "nerk", "oerk", "perk" etc. l25l . 

In the case of the 4th order gravities, statefinders may become powerful tools to constrain such theories observation- 
ally, since they enter observational relations in the higher orders of redshift z (see G71 for non-brane case diagnosis). 

Apparently, a blow-up of statefinders may also be linked to an emergence of exotic singularities in the universe 11261 . 



6. CONCLUSIONS 

We conclude the following: 

• The formulation of the fourth-order gravity theories on the brane is non-trivial because of the powers of delta 
function ambiguities. 

• Two methods were applied: 

A. Smoothing out the continuity conditions for the metric tensor at the brane; 

B. Reduction to an equivalent 2nd order theory. 

In both cases the Israel junction conditions have been obtained and they are also mutually equivalent. 

• The method of the GH boundary terms was also applied and the most general junction conditions (with no 
continuity of the scalaron and tensoron on the brane assumed) were obtained that way, too. 

• Higher-order brane gravities contain higher-order derivatives of the geometric quantities (in a Friedmann model 
it is just the scale factor) which may manifest themselves in the higher-order characteristics of expansion such 
as statefinders (jerk, snap, lerk/crack, merk/pop). 

• A blow-up of statefinders may be linked to an emergence of exotic singularities in the universe. 
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